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SDMMAEY 


The iterative interference laethod given in NACA TN No. 1252 
is applied to the solution of the following three problems concserning 
the design, of cascades: - 

(1) Determination of the shape and setting of thin (zero 
thickness) blades with a prescribed type of vortex distribution and 
total vortex strength in a cascade of given solidity for a given 
direction of the mean flow. 

( 2 ) Determination of the shape and setting of tliln blades with 
certain prescribed types of pressure distributions over ope surface 
and prescribed total, vortex strength in a cascade of given solidity 
for a given direction of the mean flow. 

( 3 ) Deterniination of the blade setting for a cascade of given 
airfoil shape and solidity that^ for a given direction of the incomiixg 
(u.pstream.) flow, will provide the front stagnation point exactly at 
the leading ed^ . 

A modification of the basic procedinre of TN No. 1252 is also 
described, in which the direction of the incoming flow, ratlier than 
the direction of tiie mean flow, is specified. 


INTROnJCTION 


In reference 1, an iterative Interference method was described 
for calculating the potential flow on an airfoil in cascade. The 
method, which malies use of charts originally employed by Betz in a 
similar study (reference 2 ), evaluates- the flow at each airfoil as 
the sum of two cemponents - that due to the uniform mean, or "free 
stream" flow, and the interference flow induced by the presence of 
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all the other airfoils of the cascade. As was indicated in 
reference 1^ such an. approach provides conslderahle flexibility and 
penults the solution, with reasonable facility, of certain cascade 
problems that would be very difficult by the usual methods that 
seek directJy^ the confonoal tranafonnation of the cascade to a 
circle! In the present paper the solutions of three such problems 
are described and an example of each is given. Two of the problems 
concern the design of thin (zero thickness) airfoil cascades having 
prescribed types of vortlolty distribution or pressure distribution 
along the blade. The other pixjblem is the detemiination of the 
blade setting, for a cascade of given airfoil shape and solidity, 
that, for a given direction of tlie inccanlng (upstream) flow, will 
provide the front stagnation point exac tly at the leading edge . A 
modification of— the basic procedure of reference 1 is also described 
in which the direction of the incoming flow, rather than the direction 
of the mean flow, is specified, 

Basic concepts and techniques are given in reference 1 , and 
detailed discussions are accordingly given herein only for those 
parts of the procedure that are not contained in reference 1. 


SYMBOIS 


Y flow velocity 

u component of flow velocity parallel to stagger line 

V ccsaponent of mean flow velocity normal to stagger line 

p blade angle, angle between chord and normal to cascade axis 
a angle between chord and velocity indicated by subscript 
stream function 
$ velocity potential 

r circulation required to provide stagnation point at trailing 
edge 

r* olroulation required to provide stagnation point at leeuilng 
edge 

7 circulation per unit etro length 
x,y coordinates of point on arc 
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0 angular position of point on aro 

p radius of circular arc 

r chord distance "between points on arc 

B arc length 

li. increment of arc length 

ajh points on arc 

O) angle of chord between points on arc 

n displeiGement normal to arc 

S angle of deviation 

E radius of tranafoamied circle 

e difference between circle and near circle angles 



Subscripts 

0 mean flow 

1 incoming flow 

1 interference on central blade due to presence of external blades 
s self— induced flow 

N ‘ normal to cascade axis 
P parallel to cascade axis 
T total 
n nose 
t tralliiig edge 
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zei"Oth, firsts second, and third approximations 
u, I upper and lower surfaces 

a,h points on arc - 

PEOBLEMS IN THE DESIGN OF TEEN BLADES 
Prohlem I 


The problem most readily solved Is the determination of the 
shape and setting of thin blades with a prescribed type of vortlcity 
distribution and total vortex strength in a cascade of given solidity 
for a given direction of the mean flow. The procedure is as follows i 
a reasonable blade shape and setting Eire assumed- for the zeroth 
approximation of the iteration process. The stream function and 
velocity potential on a particular blade of the cascade, which will 
be I’ef erred to as the central blade, are caicvilated as the sum of 
those due to the following three flows; 

(1) The given mean flow 

(2) The cascade interference flow, due to the vorticlty distri- 
bution on all the other blades 

(3) The self— Induced flow due to the vorticity distribution along 
the central blade itself 

The assumed shape is then rotated and distorted to make it 
approximate a streamline in this flow (that is, to make it coincide 
with a line along which the stream fiuiction is constant) and this 
new shape is used for the cascade blades of the next approximation. 
Through iteration of this procedure iintil distortions are too small 
to affect the flow, the desired blade shape and setting are obtained. 

Mean flow .- For convenience, the component of the mean flow 
velocity normal to the stagger line v vill be considered as unity. 
The velocity Vq of the mean flow is then known from its direction. 
The stream function to and velocity potential q a point 

x,y on the central blade are then seen from figure 1 to be; 
t 0 (Xjy) * - "Vq (x sin oq + y cos oq) j 


( 1 ) 
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where oq is the angle of attach of the mean flow with respect to 
the chord. 

Cascade interference flow ,— Tte presorihed vorticity distri- 
bution on. each blade of the cascade is replaced by a number of 
discrete vortices, and the induced stream function ilr and velocity 

potential $ are found by using the charts and methods discussed 
in reference 1. 

Self— induced flow. — The stream function induced at a point b 
of the central blade by the vorticity distribution on that blade is 
given by the integral along the blade of the imaginary part of the 
complez flow f\mction of the distribution, namely 


ts(sb) 



pst 

I r(3a) loSe^ab «^a 
Jsn 


(2) 


where y is circulation per unit arc length, and rQ_-p is the 
length of the chord between b and the variable point a. 

As the variable point a approaches the point b, the Integral, 
becomes improper. If the segaent from s-^j — n to is made 

sufficiently small, that portion of the blade can be considered a 
straight line and its vorticity uniform. Then r^b becomes 

js)j — Sa| and 7 {s) becomes a constant /(s'^j). The stream function 

induced by that segnent then is 


7(Sb) 

2« 



Us-b - 





7(sb) 





( 3 ) 


The total self— induced stream function at the point b is 


then 
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2« 


i^t 

/ 7(Ba)l°Serab 

Ja^+U2 


7(s-b) 

•‘5«~ 


lij_ - H2^10gQHj_ + ti2 - U2l°Sel^2 


(h) 


Since the function loggrat is not 3mown amlyticalJ^ for an 
arbitrarily s2xaped blade^ the integrals of equation ( 4 ) are evaluated 
numerically. 

The corresponding integration for the velocity potential requires 
special care, inasmuch as the contribution of each vortex 
element 7 (s)ds is multivalued. For the present case uniqueness may 
be provided by using the section of the blade that lies to the right 
of the vortex element as the branch cut. In which case the potentials 
contributed by a vortex element situated at a point a are defined 
by the angles shown in figure 2« For a point b to the left of 
a, the angles defixiing the potentials on the upper and lower sides, 
b^ and b^, ai’e the same. For a point b on the right side of 

a, however, it is necessaiy to differentiate between b^ and b^j 

thus the angle representing the potential on the upper side is 
designated a^abu* representing the potential on the lower 

side is designated where, from figirre 2, « 2 tt + t^ab^* 

The velocity potentials due to the entire blade at the upper and 
lower sides of a point b are then given by 

7(sa)=«abu ^a + A 7(Ba)«%hu dSg 

1 r^b 1 

= St / 7Csa>»ab; ^ 7(sa)^b^ 

■^®n Js^ 

,^b ^ rt 

/ 7(saKba ^ 7(sa)(2jr + 

Jb^ ^®b 


1 


( 5 ) 
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The average self -induced velocity potential at the point b is then 


® s(®b) 


^u(sh) 


r 


1 r* 


’^^^aKbu ^®a + ST j r(8a)<^a-bu 


St 




2 


7.(8a)<^a 


I St) 


(6) 


For convenience, a new angle £%{j is ^acw defined which is always 
measured from right to left, irrespective of the position of a 
with respect to b, and for which the following relations then hold: 


"“abu 

II 

for 

®b > ®a 

^’ab^ 

= “ab “* ^ 

for 

8b ®a J 


(7) 


The average self— induced velocity potential at a point b is there- 
fore given by substituting relations (7) into equation (6), thus 


^ sC^b) 



rCsahab <^£ 


( 8 ) 


Since the angle is not Imown ane lytioally for an arbitrarily 

shaped blade, the Integral of equation (8) is evaluated numerically. 

The choice of a circular arc shape for the initial approximar- 
tion facilitates the calculation of the self— Induced flow, especially 
when xartform vorticity is specified. The self— induced velocity 
potential and stream function for constant vorticity on a circular 
arc are. derived in the appendix. 


Rotation and distortion of the blade shape ,— The sTan of the 
stream, functions ="^0 + ''!^s ■*" 'J^i assumed blade will in 

general not be the same at every point; that is, the blade will 
not be a streamline in the ocaaplete flow. The flow crosses the 
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Talad© at each, point at an angle given by the rsufcio of the local 
noimal velocity to the local "velooity. Thtug 


6(b) ^ 

ci^r^/ba 6^y/da 


C9) 


where fi(s) is the deviation angle measured clochwise^ and n Is 

the coordinate normal to the arc at any point. In order to make the 
blade follow the streamline, the blade npaet be distorted so that 
the direction of each element is changed by this deviation angle. 

For emadl deviations, the distortion is effected by a nornnal displace— 
35iient given by the integral of the deviation angle along the blade 


As an inteimedlato step_the given shape may be rotated by the 
average angle of deviation 5 before the blade shape is distorted. 
New mean flow and interf-trenoe velocity potentials and stream 
functions are then foun'j. (the self— Induced flow remains the same, 
however) and the distortions are calculated as Just described. In 
order to minijnlze the distortions, displaoanente are taken relative 
to the displacement at the point where the calculated displacement is 
the average of the extreme displacements. 

If a finite vortlclty is required near the tips, infinite slopes 
appear at the tips (reference P« 10). The finite elopes that 
are defined by the arbitrary procedure of reference 3 oan, however, 
be adapted to the present problem. The blade shape at the tips 
given in inference 3 s- similar vortlclty distribution can be used, 
if the ordinates are multiplied by the lift coefficient of the blade 
baaed on the average tangential velocity at the tip. 

Example I .— As an example of the design procedure outlined, a 
blade was designed for a cascade of solidity 1,5 such that, with a 
mean flow direction making 6in angle of h0° 5^* with the normal, to 
the stagger line, it would have uniform vortlclty along the blade and 
a total vortex strength of 1.7321 per blade (based on vinlt velocity 
noonaal to the stagger line and unit cascade spacing). These condi- 
tions correspond to a flow coming in at 6o° to the normal and 
leaving nozmal to the stagger line. 
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A 60° circular arc (indicated "by 0 in fig. 3) vlth a chord of 
1.5 was chosen as the zeroth approximation. The blade angle was 
taken as 30 ° to the normal so that the tips were tangent to the 
Incoming and outgoing flow directions. The angle of attack of the 
mean flow is then 10°5^** Ey following the indicated procedure the 
average angle of deviation was' found to be 7 ^ 31 * . Eotating the 
circular arc through this angle and repeating the procedure described 
gave the shape designated as I in figure 3* The rotation angles for 
the second and third appr-' cimations were 0°49* and 0°9*, respectively. 
The shapes obtained are designated as II and III in figure 3 . The 
angles of deviation 6 found after rotation are plotted in figure 4 . 
The final blade angle was 36°35*« The squares of the velocities on 
the upper and lower sides are .plotted iu figure 5 , 


. Problem II 

The procedure of problem I can easily be adapted to the design 
of thin blades having cei’tain prescribed types of velocity or 
pressure distributions instead of vorticity distributions. In the 
second problem the solidity of the' cascade, the mean flow direction, 
and the total vorticity per blade are given, and a •chin blade having 
the. prescribed type of velocity distribution over one side is 
Bou^t, The procedure is somewhat similar to that of the first 
problem, A reasonable blade shape, blade angle, and vorticity 
^stribution are assumed for the zeroth approximation, ’ The average 
total velocity potential is found as in problem ! and the average 

tangential velocity is found frcaa the slope of the curve of 
plotted against the dis-tance s along the blade. This average 
tangential velocity is now plotted against s and the prescribed 
type of velocity distx’ibution on one surface is then plotted on 
the same graph so that the area between the ■two curves is equal to 
half the desired total vorticity, (The possibility of uniquely 
performing this last step determines whether the type of velocity 
distribution specified in the pTOblem is one for which a solution can 
be found by this procedure). The velocity distribution on the 
other surface may now be p_otted such -that the average— velocity 
curve falls midway between it and the curve of the velocity on the 
first surface. The area between the velocity curves for the upper 
and lower surfaces then represents the total vor'tei streng-bh, and 
the diffe^rence in ordinates at each value of s. represents the local 
vorticity on the blade. 

The vortex distribution thus detemlned is used, to find a 
first approximation -bo the shape ^d blade angle by carryiiig out 
one step of the px^scedure of p^blem I, By assuming the vortex, 
distribution to be unchanged, -the total average velocity potential $ 
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is then foimd for this new hlade and the average tangent iaJL 
velocity is detemined, A now vortlcity distribution is then found 
by the same procedure as before. The process is continued until 
further changes become Inappreciable* 

Example II .— By following the procedure outlined, a cascade of 
solidity 1.5 was designed so that the totel vortex strength per blade 
was 1.7321 in a mean flow making an angle of with the noimal 

to the stagger line, just as in the first example, but with the 
velocity on the upper surface uniform over the forward 60 percent of 
the arc and then decreasing linearly to the mean velocity at the 
trailing edge. 

The blade shape, blade angle, and vortex distribution obtained 
in example I were chosen for the zeroth approximation. The new vortex 
distribution was found by plotting the average tangential velocity, 
as in figure 6, and then finding an upper surface velocity distribution 
of the prescribed type such that the area between the two curves 
was 0.866. The vortlcity distribution was then twice the difference 
between the two ourvee. 

The average angle of deviation was found to be —1^32* and the 
shape after distortion was that designated ae I in figui'^ J. The 
rotation angles of the second and third approximation wore — 1°1* 
and 0°19’, respectively, and the shapes obtained are designated as 
II and III in figure 7« The angles of deviation & found after 
rotation are plotted in figure 8. The final blade angle was 36°21', 

The squares of the velocities on the upper and lower sides are 
plotted in figure 9. 


PEOBIU© IHVOLTING A SPECIECED INCOMIWG-ELOW DIRECTION 

Problem HI 


In reference 1 it was shown how, after a solution had been fouM 
for a given cascade in a particular mean flow, the conformal trans- 
formation of the cascade to a circle could be found; whence the 
solution for any other specified mean flow, incoming flow, or outgoing 
flow can be obtained. The present section will discuss the procedure 
for getting the solution of a given cascade directly when the incoming 
flow direction, rather than the mean flow direction, is specified. 

This problem is considered of interest because in experimental cascade 
studies the incoming flow direction is normally used as a basic 
parameter rather than the mean flow direction; furthermore, the 
discussion of this probJnm wild provide a convenient beisls for the 
discussion of the succeeding problem. 
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Mean flov atnd total vortex etrenath .-- In. the first step, as in 
reference 1, the oonfomal transformation of the isolated blade to 
a circle is determined by the methods of reference U, The flow 
field at the isolated blade is now considered as beliog oompoBed of 
three superin^josed flow fields 

(l) A uniform flow of unit velocity, normal to the stagger line ' 
(making an angle cjj with the chord) plus vortices on the blade of 

total strength Pjj which maintain tho blade a streamline in this 
flow. 


(2) A uniform flow of velocity Uq to be determined, parallel 
to the stagger line (making an angle op with the chord) plus ' _ 
vortices of strength Pp which maintain- tho blade a streamline in 
this flow. 

(3) The interference flow due to the' vortices that represent 
all the other blades of the cascade, plus vortices of strength P ^ 

which maintain the blade a streamline in this flow. The total 
vortloity on the blade will then be P,p xspjj + Pp + Pj_. 

By equation (35) of reference 4- 


Pp = 4-nRuQ sin(ap + (H) 

where E is the radius of the transformed circle and is the 

value of the difference between the circle and near-circle angles at 
the trailing edge. Similarly- ^ , 

Pj^ = 4nE sin (aj{ + £-t) = — 4i!E cos (op + e.^.) (12) 

where the relations between the angles ap • and ctjj and the magnitudes 
of the normal and parallel veloci-ty components are shown in figure 10. 

For tho third ccsnponsnt a reasonable distribution of vorticity 
along the ex-temal blades is chosen -wi-th 'a- total vortex stren^h that 
has temporarily been made unity; .ai^d the corresponding change of vortex 
strength on tlie central blade g is found by the method of reference 1. 
The ac-tual vor-teix etrong-fch induced by the external blades is the 

product of g and, the -total vortex streng-tli on eaoh blade "F^. 

The following eq,uatlon then yields the -^tal vortex .streng-fch '•within 
each approximation: 
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Fiji »• ripg'+ 4nE sin (cxp + - coe (op + (13) 

A second relation between Pip and Uq Is given by elementaiy 
cascade theory 


\lQ ■ ui - ■ 


(IM 


Simultaneous solution of eqixations (13) and (l4) then gives the 
zeroth approximation value of the total circulation and of the 
parallel component of the mean flow uq 


p 4irH sin (o£p + — 006 (op + €-t)'J 

T “ 1 _ g + 2 ttE sin (ap + €^) 

2irE [up sin (otp +• €^} — 00 s (ctj> + Cj. 
^ 1 1 “• g + 2jtE sin (<xp + 



( 1 ^) 


Vortloity distribution *- When the mean flow has been found, the 
total velocity potential on the central blade and the vorticlty dls~ 
tribution are found as in reference 1. With the uao of this vorticlty 
distribution to calculate g, the entire process is repeated to 
determine a new Pp and vorticlty distribution. The procedure is 

continued until further changes are inappreciable. The velocity 
distribution on the bJade is then found by adding the velocity due to 
the interference to that due to the mean flow, as is done in 
reference 1, or it may be found directly by differentiating the 
potential with respect to the distance along the surface. 

Example III .— The potential flow was foxond for a cascade of 
solidity 1.5 (given shape and blade angle) with incoming flow at 
45° to the stagger line. The blade shape ires that derived in 
example I and the blade angle was 21^35’ which gives the same angle 
of attack with respect to the incoming flow as in example I, 

Uniform vortex distribution was assumed for the zeroth approxi- 
mation, , The total vortex strength Pp found for that distribution 

was 1,2302 and the mean flow was at 21°3*. Successive approximations 
gave the following values; 
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Approximation 

i 

Vortex strength 

Mean flow direction 

0 

1 

i 

1.2302 

21° 3* 

I 

1.2530 

19° 40* 

II 

1.2457 

20° 40 ‘ 

III 

1.2472 

20° 38* 


The squares of the velocities on the upper and lower surfaces 
are plotted in figure 11. 


Problem XV 

The procedure of probl«a III can be readily extended to determine 
the stagger angle, for given incoming flow direction, at which an 
airfoil in cascade is at the "ideal” angle of attack. 

Ideal angle of attack condition .— The ideal condition for an 
airfoil is one for which there is a steignation point at the nose, or 
for zero thickness, the condition for which air enters tangentially 
at the leading edge. With regard to the flow in the plane of the 
circle to which the airfoil t3:anBforms, it is the condition for 
which the same vortex ceoicels the velocity at both the leading-edge 
and trailing-edge points. 

The strength of the vortex at the center of the circle 

which cancels the velocity at the leading-edge point is, in analogy to 
equation (13), 


r'l =r ‘jS* + 


— Uq sin (ttp + + cos (op + 



(16) 


The procedure of reference 1, modified to cauacel the induced velocity 
at the leading-edge point, is again used to find the factor for this 
vortex strength induced by the external blades g*. 

The ideal angle of attack condition is then the condition at 
which rp = r*Qj or Pq, — «s 0, where the total induced vortex 

strengths P^ and P'rp are found by equations (13) and (16). 
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Determination of blade angle .— A Made angle Pj Is assumed 

for the first step and the procedure of the first approximation of 
example III ie carried out to obtain and and the vortlcity 

distribution corresponding to one of these two values of the total 
vortex strength Is found. With this vortlcity distribution assumed 
on the external airfoils, the calculation is repeated for a second 
blade angle Pjj, which ie choson greater or less thsui Pj accordingly 

as Pjip is greater or less than P*j,j,‘ interpolating between 

or extrapolating from the results for those two .values of P, a 
third value of P is found for which P^ — P*rj, should be very 
nearly zero. A calculation at this value of P either should verify 
that P^ — P’y is practically zero or should provide the data for 

a more accurate interpolation or extrapolation. 

In the procedure as Just described, only one approximation is 
made for each blade angle; that Is, for each p the vortlcity 
distribution found for the preceding p (using either P^ or P 

is usod for the external airfoils. This method should, in general, 
suffice for satisfactory convergence; in any case, no more than two 
approximations for each angle should be required. 

Example IV .— Tbe blade angle was found for the blade derived in 
example I'such that it would be at -the ideal ang;le of attack in a 
oeiacade of solidity 1,5, with incoming flow to the stagger line. 

The initial blade angle and vortex distributions were those found 
in example III. Thus Pj = 21°35* and rjip = 1.2h72. By 
equation (I 6 ), r*ii' = 1.3923, Therefore r*jy = —0,lh31, 

A second calculation, with Pjj = 20%5* gave T j — P * j jip * 0,0130, 
Interpolation between the two results indicated that P^n — P’y would 
be zero at 20 °Jj 9'. A final calculation with Pjjj = 20°hp^ verified 
this fact and gave a vortex strength of 1,2491, 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va , , Janueiry l4, . 194? 
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APPENDIX 

SEIE -INDUCED FLOW FOE A CIRCULAE ARC NITH 
CONSTANT VOETICITI 


The stream function induced at point h by an element ds at 
point a is (fig;. 12) 





loBer-ab 


(Al) 


With the chord and the arc element espiressed in terms of tho angular 
position of the points, this eq.uation becomes 


lts(0b) 


2Jt' 


loge 


2p sin 



(A2) 


vhere 0 is the angular position of the points on tli© arc, and 

®b " ®a ^ 

p the radius of tho arc. Or, with 


i^Q(0b) 



(A3) 


Eq.uation (A3) is integrated witii respect to from. to 9^^ 

and tho first three terms of the series expansion of the integral 
are retained. The result of tliis integration, after all terms which 
do not contain 6 -^ have been omitted (since they add only constants 
to the stream function) is 


^s(®b) 




loSe (®t ” ®b) 



lo^ (®b ■ ®n) 




(e 


n - ©t) + ©b (®t^ - ®n^)l ^ 


(A4a) 
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By adding oonstante necessary to conplote the 6q,uare in the last 
term, this equation may ha put in the more convenient form 


= - f •! 


^9t - Qh> 


lose (0t - ©h) + 





r-* — 

1 /' N 

V wy 



(®h " ®n) 


(A4b) 


In order to find the calf -iiiduced average velocity potential 
®g for the circular arc, it is necessary to find the expression for 
tlie angle <%h previously defined. From figure 12 it is seen 
tiiat 


«%h 


« 

2 




(.A5) 


Substituting this expression into equatioj; (■?•) and di’opping.oAditivo 
constants gives the desired velocity potential 
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Fig. 1 



Figure 1.- Potential of the mean flow at a point on the blade. 


Fig. 2 
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Figure 2. - Armies defining seK-induced potentials on upper and lower 
surfaces at point b due to vortex at point a. 





NACA TN No. 1254 


Fig. 3 



Figure 3.- Shape aiKi settii^ of blade used for zeroth approximation of 
example I, and of blades derived in the subsequent approximations. 
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NACA TN No. 1254 


/z/rsf approx /mahon 
)G opprox/mpt 



concf approximation 
' opprox/mafion 



0 2.0 40 60 60 /OO 


Disionce along bbde surface, pence nf 

NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 


Figure 4.- Deviation angles between arc and streamline for the three 
approximations of example I, showing rate of convergence. 
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Figure 5.- Pressure distribution on airfoil derived in example I. 




Fig. 6 
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Figure 6. - Average and top surface velocities for example II. 
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Fig. 7 



Figure 7.- Shape and setting of blade used for zeroth approximation of 
example n, and of blades derived in the subsequent approximations. 
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Figure 8.- Deviation angles between arc and' streamline for the three 
approximations of example n, showing rate of convergence. 
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Fig. 9 
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Figure 9. - Pressure distribution on airfoil derived in example n. 
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Fig. 11 
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Figure 11.- Pressure distribution on airfoil in example HE. 



Fig. 12 


NACA TN No. 1254 


r 



NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 


Figure 12.- Definitions of angles and distances for derivation of 
self -induced potential and stream function on circular arc. 


